Control of Wannier orbitals for generating entanglement of 
ultracold atoms in an optical lattice 
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We propose a method for controUing the quantum state of ultracold atoms in an optical lattice 
by using current atom manipulation techniques, which can generate high-fidelity entanglement re- 
sources. We utilize higher Wannier orbitals as a controllable and accessible environment, and then 
realize a tunable Ising interaction between atoms in the lowest orbital to generate multipartite en- 
tangled cluster states. We can enhance the fidelity by employing a post-selection scheme based on 
spectroscopic measurements of states in the environment. Precise numerical simulations using real- 
istic parameters suggest that our method offers significant advantages for high-fidelity entanglement 
generation toward scalable measurement-based quantum computation. 

PACS numbers: 03.67.Lx, 03.67.Bg, 37.10.Jk, 71.10.Fd 
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Measurement-based quantum computation (MBQC) 
[J, 0] is a prominent method for scalable quantum in- 
formation processing because we can perform a universal 
quantum computation with only successive single-qubit 
measurements after preparing large entangled resources, 
such as multipartite entangled cluster states. Ultracold 
atoms with pseudospin {e.g. hyperfine) states in an op- 
tical lattice are promising candidates with which to im- 
plement MBQC [sj. This clean and isolated system will 
sustain quantum coherence for a long time [3|. State- 
of-the-art experiments using coupling between lights and 
atoms demonstrate that a single atom can be detected 
and manipulated by addressing the lattice sites individu- 
ally d, Q , suggesting substantial progress on single-qubit 
measurement techniques. The other key issue, the mul- 
tipartite entanglement generation, has been discussed in 
theoretical studies and has also been demonstrated 
experimentally fid\. However, high fidelity and scalabil- 
ity, which are important features in the implementation 
of MBQC, have not been fully confirmed even theoreti- 
cally. Especially, change of the atomic orbitals and mul- 
tiple site-occupancies, which have not been sufficiently 
taken into account, would indeed cause non-negligible er- 
rors in scalable quantum computing. 

In this paper, we propose to exploit Wannier orbitals as 
controllable quantum states for generating high-fidelity 
entangled cluster states of ultracold fermionic atoms in 
an optical lattice. To make this control precise, we treat 
Wannier orbitals in the ab initio manner. In our method, 
atoms in the lowest orbital are chosen as qubits, and the 
extra-Hilbert space that originates from higher orbitals 
serves as a controllable and accessible environment. By 
controlling the coupling between the qubits and the en- 
vironment, we can design interactions among the qubits. 
Specifically, we create a tunable Ising interaction in ad- 
jacent sites to generate entangled cluster states. The 
fidelity can be enhanced by performing measurements on 



states of the environment followed by post-selection de- 
pending on the resulting outcomes. Moreover, substan- 
tial advantages as regards scalability can be obtained by 
our pair-wise entanglement generation scheme. Precise 
numerical simulations involving an exact diagonalization 
confirm that a combination of the above tricks allows us 
to generate very high-fidelity entanglement with current 
experimental technologies. The present method is appli- 
cable to generating one, two, and three dimensional (ID, 
2D, and 3D) cluster states, and thus is suitable for fault 
tolerant MBQC schemes 

We begin by explaining our system. We consider 
fermionic atoms with two spins in an optical lattice. The 
non-interacting Hamiltonian of the atoms is given by 
Ho = ~Y.-^^-^-hVo[cos{x.Ja)+coa{yi/a) + cos{zi/a)], 
where Vq is the lattice depth, a is the lattice dis- 
tance, and M is the mass of the atoms. We solve 
Ho numerically, and calculate the eigenvectors and their 
Fourier transformation, namely, Bloch and Wannier or- 
bitals. These Wannier orbitals are used to derive the 
full Hamiltonian, H = Ho -\- Hint, hi the second quan- 
tization framework. Here, Hint is the contact interac- 



tion term given by 
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J2i j - ^j), where as is 



the scattering length. The second quantized expression 
of H is the multiorbital Hubbard Hamiltonian written 

&S H ~ '^{ij) EcTQ '^aaC^^a'^jcra + Ei Ecrct ^a'^iaa'^i'^a + 

where cl^^{ciaa) are the creation (annihilation) opera- 
tor of a fermion with a spin o'(=^,4) and an orbital 
a{= 1, 2, • • • ) at the ith site, and the subscript (ij) means 
the summation over the nearest-neighbor sites. Parame- 
ters, on-site potentials hopping matrices to adjacent 
sites Ja/3 , and interaction strengths Uap , are determined 
by ab initio calculations. We note that inter-orbital hop- 
ping matrices vanish due to orthogonality of the set of 
the Bloch functions, i.e., Ja^js = 0. We assume that the 
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FIG. 1: (a) A schematic diagram of a lattice potential modu- 
lation induced by the additional two-site period potential, (b) 
and (c) The parameters of the additional potential Hamilto- 
nian Ti' for Vo = 15i5r and Vq ~ 4Er. Onsite potential terms 
e'i^a are real number, Ime^ „ = 0. All elements of J^^ except 
for ImJ{2 are close to zero. 
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FIG. 2: (a) and (b) Second order virtual hopping processes 
for two atoms with different and same spins, respectively. 
Blue solid lines represent the energy levels of the lowest and 
the second lowest orbitals. The gray dotted line shows the en- 
ergy shift caused by the interaction, (c) Resonant condition of 

the Ising interaction with a strength Jiaing oc ^^^^^ — ^ . 



trapping potential is absent from our model as recently 
demonstrated in the experiment |15i] . 

In what follows, we define single atoms with a =t (i) 
occupying the lowest orbital a = 1 at the ith site 
as the qubits, and express its state by |0)i(|l)i). As 
we will see later, the other states in the higher or- 
bitals or with multiple occupancies serve as the control- 
lable and accessible environment that plays an essential 
role in our method. It is useful to mention here the 
general strategy for creating the cluster state |CS) f^]: 
first, initialize qubits as the product of superposition 
states = (|0)i + \l)i)/\/2, which are efficiently pre- 
pared using the band insulating transition of fermions 
[la ]: then, apply a unitary transformation described by 
the time evolution with an Ising interaction, namely, 

ICS) = eM-Kr/h)Jl,insE^^!^^+l}Yl^\+)^ ^Ith the 

characteristic time r = 47rfi,/ Jising- We will show that 
the unitary transformation can be implemented by sim- 
ply modulating the optical lattice potentials. 

Now, we explain our implementation scheme for re- 
alizing the Ising interaction. In addition to the origi- 
nal lattice potential, we apply an extra lattice potential 
with a two-site period given by Vq cos(a;/2a -|- 9), where 9 
is the relative phase between the original and the ad- 
ditional potentials. Here, we chose 9 = tt/2 and set 
it along the a;-direction. Figure [T](a) shows schemati- 
cally how the combination of these two potentials mod- 
ifies the lattice landscape. The Hamiltonian originat- 
ing from this additional potential is written as "H' — 

We evaluate 



the parameters e'^ ^ and with the same ab initio cal- 
culations as mentioned above. In Fig.[T](b), the on-site 
potential e exhibits staggered modulations as expected 
from the total potential shape shown in Fig[T](a). On the 
other hand, as regards hopping matrices J'^ we find an 
interesting and useful feature as shown in Fig.[l](c). Only 
the inter-orbital J[ ^2 has a large magnitude, while the 
other elements J[ are nearly equal to zero. This is be- 
cause the orthogonality of the Bloch functions is broken 
by superposing the two-site period potential. In addi- 
tion, this modulation is commensurate with the original 



lattice so that strong inter-orbital coupling is introduced 
between orbital- 1 and 2. Such a striking situation is re- 
alized when 9 — 7r/2 -I- ^tt with integer ^. As 9 deviates 
from these values, some elements of become non- 

negligible. 

Next, we roughly explain the physical mechanism by 
which the two-site period potential leads to the Ising in- 
teraction, based on the second order perturbation theory 
with specific initial states such as Mott insulators. As 
is well known, such perturbative processes originating 
from Jii induce Heisenberg type spin interactions with a 
strength JHciscnbcrg oc Jli/Uii. We now consider similar 
processes in terms of the inter-orbital hopping J'12 with 
a large amplitude. Note that the z-dependence of J[ 12 
is currently neglected for simplicity. Figure [2](a) and (b) 
schematically represents the perturbative processes for 
two localized atoms with different and same spins, re- 
spectively. The top panels correspond to the initial Mott 
states. The middle panels show the virtual intermediate 
states. It is clear that, for the atoms with different spins, 
the inter-orbital interaction U12 causes the energy shift 
shown by the gray dotted line, while there is no such en- 
ergy shift for the atoms with same spins due to the Pauli 
principle. The bottom panels are the final states, which 
are exactly the same as the initial states. This is a conse- 
quence of an energy conservation law that prohibits spin- 
exchanging processes with a finite energy difference due 
to the Bloch bandgap A — £2 — £i + £'i 2~^i+i i- Thus, we 
can conclude that the additional potential results in the 
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Ising interaction with a strength Jising oc X+U12 ~ a 
via the perturbation processes. 

Here we remark on the controllability of our Ising in- 
teraction. We can switch off the interaction at our will 
by removing the additional lattice potential. More im- 
portantly, the interaction strength is increased if the res- 
onant condition A -I- \J\2 ~ is satisfied as shown in 
Fig.[2](c). We can achieve this resonant condition by con- 
trolling J7i2(oc as) via a Feshbach resonance 17| and also 
by tuning £^ ^ via the change in VJ,'. The resonant con- 
dition of the Ising interaction allows us to reduce the 
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FIG. 3: Simulation results for the 2-qubit system with Vo = 
lOEr and Vq = 6.2Er. (a) Fidelity F as a function of time 
r. (b) The number of doubly occupied sites D and the atoms 
with the second lowest orbital Ai'2nd- The two curves over- 
lap each other, (c) Fidelity with the post selection Fps and 
success probability Psuc 
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FIG. 4: (a) Overall diagram of a linear potential modulation 
in addition to the lattice modulation shown in Fig.[T](a). The 
dotted line provides a stepwise structure to guide the eyes, 
(b) A schematic view of the pair-wise scheme, where the Ising 
interaction is resonantly enhanced only inside the unit cells. 



time required for preparing cluster states. Note that the 
large magnetic interaction is usually difficult to realize 
via the perturbative mechanism without such a resonance 
scheme. 

The above discussions should be carefully verified by 
calculating the actual dynamics of atoms without any 
perturbative approximations. We therefore numerically 
simulate the time evolution of the present system written 
as \4>{t)) = exp{i{T/h){'H + V.')} Yii using an exact 
diagonalization method. It is found that certain effects 
not included in the simple second-order perturbation the- 
ory degrade both fidelity and scalability of our method. 
To overcome these problems, we adopt two more strate- 
gies in addition to the resonance scheme, namely, the 
post-selection scheme and the pair-wise scheme. In the 
following, we will explain in detail how each scheme of 
our proposal works well on the basis of numerical simu- 
lations. 

First of all, we examine our resonance scheme by sim- 
ulating a 2-qubit system of ""^K atoms with the follow- 
ing realistic parameters; a = 413 nm, as = —50 nm, 
Vq = lOEr and Vq = 6.2Er, where i?,- is the recoil en- 
ergy. Note that the third lowest and the higher orbitals 
are negligible because energy levels of them are far off- 
resonant from that of the lowest orbital. We calculate 
the fidelity F = |(CS|(/>(t))|2, where |CS) is equivalent 
to the Bell pair in 2-qubit entangled states. In Fig.[3](a), 
the calculated fidelity F is plotted as a function of time 
T. A high-fidelity cluster state of F = 0.997 is obtained 
at around 1.5 ms. We thus confirm that our scheme can 
realize the large Ising interaction enough to generate a 
cluster state in a short time. 

Secondly, we further explore the characteristics of the 
F curve shown in Fig.[3](a). We can see two types of os- 
cillations with periods of about 3 and 0.5 ms. The longer 
period oscillation obviously corresponds to the time evo- 
lution driven by the Ising interaction. On the other hand, 
the shorter period oscillation suggests the existence of 
some additional dynamics during the entanglement gen- 
eration. To clarify this, we plot in Fig.[3](b) the number 
of doubly occupied sites D and the number of atoms in 
the second lowest orbital A^2nd as a function of time. Two 
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FIG. 5: Simulation results using the pair-wise scheme for 4- 
and 6-qubit systems with Vo = 18Er and Vq = AEr- (a) Fi- 
delity F and Fps, and success probability Psuc (b) Rescaled 
quantities P^/", Ppg", and PsJ", where the curves for differ- 
ent n — 4 and 6 overlap each other, (c) Same quantities with 
the unit cell shift for n = 4. The gray dashed line represents 
the time r = 3.7 ms at which we move the unit cells. 



quantities D and A^2nd agree exactly and the oscillation 
period of them is 0.5 ms. These results mean that the 
actual excitation of atoms occurs between the 2-qubit 
state and the intermediate state as depicted in Fig.[5](a). 
The shorter period oscillation of F curve is attributed 
to this excitation. Note that the intermediate state is 
assumed to be virtual in the simple perturbation theory. 
The above dynamics beyond this assumption are induced 
as a trade-off for the large Ising interaction in our reso- 
nance scheme. 

Here we explain that the fidelity can be enhanced by 
eliminating the intermediate state via spectroscopic mea- 
surements. Specifically, we employ the following post- 
selection scheme: first, detect D by using, e.g., photo- 

' j, or measure iV2nd with. 
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association spectroscopy 
e.g., orbital blockade spectroscopy [20|; after that, deter- 
mine the operation is a success (failure) by the measure- 
ment outcome suggesting the absence (existence) of the 
intermediate states. In Fig.[3](c), we show fidelity with 
the post selection Fps and the success probability Psuc 
obtained by the same numerical simulation method. As 
the result of the post selection, the intermediate states 
are removed, and a shorter period oscillation disappears 
from the fidelity curve. We can thus achieve an extremely 
high fidelity of Fps '-^ 1 in return for the failure proba- 
bility of 1 - Psuc ^ 0-002. 

Thirdly, we simulate n-qubit systems for n > 2 so as 
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to discuss scalability. When we apply our scheme to 4- 
qubit system in a similar way to that described above, 
we obtain the fidelity of Fps ^ 0.71 with Pguc ^ 0.931. 
This crucial decrease in the fidelity is attributed to ex- 
tra excitation processes that are not included in those 
shown in Fig.[5](a) and (b). For instance, in the inter- 
mediate state shown in Fig.[2](a), the excited atom does 
not always go back to the initial site, but instead further 
moves to the third lattice site. Such extra excitations 
may occur except in 2-qubit systems. 

To suppress the extra excitations, we divide the system 
into a set of effective 2-qubit systems, and then indepen- 
dently generate an entanglement in each 2-qubit system. 
We realize this pair-wise scheme by adding a linear po- 
tential gx with a gradient g as demonstrated in the ex- 
periment [21]. Figure |4](a) schematically represents the 
linear potential modulation. Here we find a characteris- 
tic stepwise structure. As shown in Fig.|4](b), we define 
each step consisting of two lattice sites as a unit cell, and 
properly tune the resonant condition A -I- U12 ~ sat- 
isfied only inside the unit cells. Consequently, the Ising 
interaction is made pair-wise, and thus entanglement of 
two qubits in each unit cell can be generated simultane- 
ously and separately. 

Using the pair-wise scheme mentioned above, we per- 
form numerical simulations for n > 2, and again dis- 
cuss scalability. Here, we calculate the fidelity F be- 
tween a product state of Bell pairs in unit cells and 
|(/)(t)). In Fig.[S](a), we show results for n = 4 and 6 with 
Vo = 18Er and Vq = AEr- High fidelity cluster states for 
n = 4 (6) are generated with or without the post selec- 
tion as Fps ~ 0.999 (0.998) with P^uc ^ 0.987 (0.980) or 
F ^ 0.973(0.960), respectively. With the ideal pair-wise 
scheme, all of these quantities should be scaled as, e.g., 
F = where the tilde represents quantities per unit 

cell. In Fig.[5](b), we show rescaled quantities P^/", Fpg" 

and Psu'c"- For all of these quantities, the rescaled curves 
for 71 = 4 and 6 overlap each other, suggesting that the 
Ising interaction is almost perfectly pair-wise. Here we 
have obtained the unit cell fidelity of Fps ~ 0.999 with 
Psuc ^ 0.993 or P - 0.987. These values reflect that 
our pair-wise scheme strongly suppresses the extra exci- 
tations. 

We further extend our scheme for generating multipar- 
tite entangled states over whole ID lattice sites. After 
the set of Bell pairs are prepared as mentioned above, 
we move unit cells toward one site right or left by chang- 
ing the sign g or the relative phase — 7r/2 — >■ 37r/2, 
and then again perform pair-wise entanglement genera- 
tion. In Fig.[5](c), we show simulation results for n = 4 
with the same parameters as above, where unit cells are 
shifted at time r = 3.7 ms indicated by the gray dashed 
line. By combining our several techniques, we finally ob- 
tain very high-fidelity cluster states, Pps = 0.998(~ -fps) 
with Psuc = 0.979(- Pf^J or P ~ 0.961(- P^). An 



important point of this result is that we succeed in 
suppressing the effects of uncommutative interactions, 
such as the Heisenberg interaction. Namely, we can set 
-'Rising > -'^Heisenberg in our scheme. 

Finally, we discuss the generation of 2D and 3D cluster 
states. So far, we have discussed the ID entanglement 
generation, where the inter-orbital hopping with a 
large magnitude is only induced along the x-direction. 
Importantly, this directivity helps us to suppress the ex- 
tra excitation processes along the y- and z-directions. 
Thanks to the directivity and commutativity of our Ising 
interaction, the successive ID entanglement generation 
by changing the direction allows us to generate 2D and 
3D cluster states, which will be useful for fault tolerant 
MBQC [11-13]. Note that the failure event of the post 
selection can be regarded as loss of qubits. Therefore, 
our scheme, which enhances the fidelity in return for the 
increased losses, is suitable for the fault tolerant MBQC 
that is robust against losses [l^ . 

In summary, we propose a method for controlling the 
quantum state of ultracold atoms in an optical lattice. 
Utilizing higher Wannier orbital states, we realize a tun- 
able pair- wise Ising interaction between atoms in the low- 
est orbital. This scheme can generate high-fidelity entan- 
gled cluster states, which is useful for scalable MBQC. As 
a prospect, we note that our tunable magnetic interaction 
can be employed for quantum simulations of magnetism, 
for instance, demonstration of the Neel transition. We 
also remark that our basic concept, Wannier orbital con- 
trolling, can be realized in both fcrmionic and bosonic 
systems. Finally, we expect that our scheme is not lim- 
ited to only these applications, but will also find a wide 
range of uses with respect to quantum state manipula- 
tions and quantum information processing. 

We thank Y. Takahashi, K. Azuma, Y. Matsuzaki, and 
Y. Tokura for valuable discussions. 
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